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Coffee  Break/  Pause  cafe 

P.  Orlandi  (Urdv.  La  Sapienza,  Roma):  3d  simulation  of  free  vortex  rings  and 
rings  interacting  with  a  wall 
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J.P.  Chollet  (LEGI,  Grenoble):  Turbulent  eddy  structures,  combustion  and  chemical 
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autour  de  structures  cylindriques  par  une  methode  de  superposition  differences 
-particules. 

Closing  address  /  Cloture 
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(1)  Laboratory  for  Aero-  and  Hydrodynamics,  Delft  University  of  Technology,  The 
Netherlands;  (2)  Faculty  of  Applied  Mathematics,  University  of  Twente,  The  Netherlands 

THE  VORTON  METHOD: 

THEORY  AND  APPLICATIONS 

The  vorton  method  is  a  vortex  method  in  which  3-D  vortex  flows  are  dis¬ 
cretized  by  means  of  singular  vortex  structures  (delta-functions),  like  in  the  2-D 
point- vortex  method  [1],  However,  a  continuous  ’’dipole”  part  has  to  be  added 
in  order  to  make  the  vorticity  field  divergencefree. 

Departing  from  Helmholtz’s  vortex  deformation  equation,  the  dynamics  of 
the  vorton  parameters  (location  and  intensity  of  vorticit,)  is  described  by  a  set 
of  nonlinear  ordinary  differential  equations  for  vorton  displacement  and  vorton 
deformation.  A  review  of  vorton  theory  is  given  in  [2}  and  [3]. 

Several  aspects  of  vorton  theory  have  been  or  will  be  investigated:  1)  motion- 
invariants:  the  vorton  dynamics  should  show  conservation  of  variables  such  as 
total  kinetic  energy  and  total  helicity;  2)  convergence:  the  vorton  method 
solutions  should  converge  towards  the  exact  smooth  solutions; 

Some  elementary  vorton  simulations  have  already  been  performed.  If  our 
newly  derived  vorton  equations  (submitted  to  Physics  of  Fluids)  appear  to  give 
reliable  and  internally  consistent  simulations  (e.g.  of  circular  and  elliptical  vor¬ 
tex  rings),  it  could  be  applied  to  the  study  of  coherent  structures  in  a  turbulent 
boundary  layer,  the  organized  vortical  structures  in  turbulence.  For  that  pur¬ 
pose,  investigation  into  the  ’’chaotic”  dynamics  of  vortons  may  be  helpful.  A 
convenient  visualization  method  will  be  developed  for  this  purpose",  using  AVS. 
This  also  gives  us  the  possibility  to  study  important  scalars  like  helicity  density. 

Meanwhile,  it  is  tried  to  derive  vorton  equations  by  making  use  of  the  fact 
that  the  3-D  Euler  equations  have  a  so-called  Poisson-structure. 

Along  with  the  singular,  or  "hard”,  vortons,  we  also  consider  the  theory  of 
smoothed,  or  "soft”,  vortons. 

An  experiment  is  set  up  in  order  to  study  the  behaviour  of  circular  and  ellip¬ 
tical  vortex  rings.  Experimental  results  will  be  compared  to  vorton  simulations. 

[1]  M.I.  Aksman,  E.A.  Novikov,  "Reconnections  of  vortex  filaments”.  Fluid 
Dynamics  Research  3(1988). 

[2]  F.  Alkemade,  "The  world  of  vortons;  origin,  theory,  applications,  prospects”. 
Laboratory  for  Aero-  and  Hydrodynamics,  Delft  University  of  Technology,  1992. 

[3]  F.  Alkemade,  E.  van  Groesen,  "A  new  set  of  dynamical  equations  for  3-D 
vortex  particles  (vortons)".  (submitted  to  Physics  of  Fluids). 


Computational  investigation  of  vortex  shedding 


Christopher  Anderson 


Department  of  Mathematics 
University  of  California 
Los  Angeles,  CA  90024 
USA 


In  this  talk  a  high  order  finite  difference  method  for  the  solution  of  the  incompressible 
Navier-Stokes  equations  in  vorticity  stream  function  form  will  be  presented.  We  will  then  discuss 
how  the  numerical  results  obtained  with  this  method  have  been  used  to  investigate  the  inaccuracies 
associated  with  the  use  of  the  time-dependent  Prandtl  boundary  layer  equations  to  model  vortex 
shedding  -  a  procedure  implicit  in  vortex  sheet  /  vortex  blob  methods.  We  will  also  show  how 
alternate  computational  strategies  can  reduce  these  inaccuracies. 


Viscous  Splitting  of  the  Navier-Stokes  Equations  with 

Boundaries 


J.  Thomas  Beale 
Mathematics  Dept. 
Duke  University 
Durham.  NC  27706  USA 

June  12,  1992 


We  consider  the  approximation  of  the  Navier-Stokes  equations  for  incompressible  flow 
by  fractional  time  steps,  without  spatial  discretization.  Assume  there  is  a  smooth  solution 
in  a  bounded  domain  in  two  or  three  space  dimensions,  with  velocity  zero  on  the  boundary. 
We  approximate  by  alternating  partial  steps:  In  the  first,  we  solve  the  inviscid  Euler 
equations  with  normal  velocity  zero  on  the  boundary.  In  the  second,  we  solve  the  linear 
Stokes  equations  with  zero  velocity  on  the  boundary,  even  though  the  initial  state  for 
this  step  has  tangential  velocity  at.  the  boundary  created  during  the  Euler  step.  In  joint 
work  with  Claude  Greengard.  of  the  IBM  T.  J.  Watson  Research  Center,  we  have  shown 
that  this  approximation  is  first-order  accurate,  that  is,  for  small  time  step,  the  error  has 
a  bound  proportional  to  the  time  step.  The  validity  of  such  a  fractional  step  method 
or  splitting  is  presumed  in  some  computational  methods  for  viscous  flow,  including  the 
transport-difFusion  method  and  vortex  methods.  In  proving  the  result,  it  seems  necessary 
to  work  with  low  norms  because  of  the  inconsistency  during  the  Stokes  step  between  the 
initial  state  and  the  zero  boundary  condition. 


Resume  pour  le  congres  ”Ecoulements  Tourbillonnaires  et  Methodes  Numeriques”,  ayant  lieu  du  15  au 
19  juin  a  Grenoble 
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Resume 


L’objectif  de  cette  etude  est  1’analyse  des  mecanismes  d ’instability  dans  un  jet  plan  compressible. 
La  resolution  de  ce  type  d’ecoulements  compressibles,  instationnaires,  par  simulation  numerique 
directe  necessite  un  schema  numerique  capable  de  prendre  en  compte  des  variations  brutales  des 
variables  sans  trop  diffuser  ni  produire  d’oscillations  numeriques  importantes  au  voisinage  de 
ces  fronts. 

Notre  choix  s’est  porte  sur  la  methode  de  Transport  pax  Flux  Corriges  developpee  par  BORIS 
et  BOOK  (1).  C’est  une  methode  explicite,  de  difference  finie,  qui  assure  la  positivite  et  la 
monotonie  de  certaines  grandeurs  physiques.  Sa  precision  a  l’ordre  4  en  phase  la  rend  tout 
particulierement  performante  dans  les  regions  de  fortes  discontinuity  et  de  chocs.  Elle  effectue 
en  tout  point  une  ponderation  entre  une  solution  du  premier  ordre  (monotone  mais  diffusive) 
et  celle  du  deuxieme  ordre  (plus  precise  mais  dispersive),  afin  de  garantir  la  positivite  et  la 
monotonicite  de  la  masse  volumique.  Son  principe  est  d’ajuster  localement  le  coefficient  de 
diffusion  en  fonction  du  profil  de  vitesse  et  de  la  distribution  de  la  variable  convectee  pour 
conserver  a  la  fois  la  stabilite,  la  monotonicite  et  la  precision  du  schema.  C’est  une  methode 
essentiellement  non  lineaire. 

Pour  tenir  compte  de  l’influence  amont-aval  en  subsonique,  deux  types  de  conditions  aux  limites 
sont  envisages:  le  premier  travaillant  avec  les  relations  dites  ”  caxacteristiques”  sur  les  derivees 
des  variables,  le  second  travaillant  directement  sur  les  variables  caxacteristiques  ou  Invariants 
de  Riemann. 

Le  caicul  du  jet  a  ete  fait  a  un  nombre  de  Reynolds  base  sur  le  diametre  de  10000.  L ’instability 
des  structures  tourbillonnaires  organisees  a  ete  bien  captee  comme  l’illustre  les  isomachs  figure 
1  et  les  isovaleurs  de  la  composante  transversale  de  la  vitesse  figure  2.  L’analyse  de  l’ecoulement 
a  permis  de  distinguer  trois  regions:  La  premiere,  proche  de  la  sortie  ou  la  vitesse  sur  1’axe  reste 
sensiblement  constante  et  egale  a  sa  valeur  a  remission,  dite  ”Zone  a  potentiel”;  la  seconde,  loin 
a  l’aval  ou  les  proprietes  du  jet  obeissent  a  des  lois  de  similarity  relativement  simples  et  bien 
connues;  et  la  troisieme  "Zone  d’interaction”  qui  sert  de  transition  entre  les  deux.  L’analyse  spec- 
trale  realisee  sur  le  signal  temporel  de  vitesse  et  de  pression  a  permis  de  retrouver  la  frequence 
caracteristique  d’emission  des  structures  organisees  (  cf  figure  3). 
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Fig.l:  Isomachs 
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Analysis  and  numerical  simulation  of  turbulent  flows  using  wavelets 


C.  Basdevant 

LMD-ENS 
Rue  Lhomond 
75005  Paris 


The  wavelet  transform  decomposes  a  signal  into  contributions  localized  both  in  physical  space  and 
spectral  space.  For  turbulent  flows  which  exhibit  a  large  spatial  variability,  it  can,  much  better  than 
the  classical  Fourier  transform,  separate  active  regions  such  as  vortices  from  the  background 
flow.This  allows  detailed  space  scale  analysis;  with  the  wavelet  transform  we  are  able  to  define 
locally  in  space  and  energy  spectrum  and  locally  in  space  and  scale  a  Reynolds  number.  This  will 
be  illustrated  by  several  examples.  On  the  other  hand,  the  wavelet  transform  can  be  used  for 
information  compression;  indeed  most  of  the  non  vanishing  coefficients  of  the  decomposition  of  a 
field  are  concentrated  in  active  regions,  when  intermittency  is  strong,  omly  a  few  coefficients  are 
needed  to  represent  the  field.  This  property  is  the  basis  of  adaptive  numerical  methods  using 
wavelets  for  numerical  simulation  of  turbulent  flows,  we  will  review  some  of  these  methods. 


Magnetization  methods 
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There  is  a  canonical  structure  for  incompressible  fluid  flow  which  is  canonical  in  any  number  of 
space  dimensions.  We  describe  lagrangian  numerical  methods  based  on  the  canonical  hamiltonian 
equations  and  present  numerical  results  for  the  two  and  three  dimensional  cases.  In  order  to 
describe  this  canonical  hamiltonian  structure  a  new  variable  is  introduced  whch  we  call  the  velocity 
magnetization.  The  velocity  magnetization  has  units  of  velocity,  but  is  related  to  the  vorticity  of  the 
fluid.  A  derivation  of  the  continuous  equation  of  motion  for  the  velocity  magnetization  in  three 
dimensions  is  presented.  The  continuous  equation  of  motion  is  discretized  using  a  regularization 
kernel  in  such  a  way  that  the  hamiltonian  structure  is  preserved  in  the  discrete  system  and,  more 
importantly,  the  discrete  dipoles  obtained  do  not  have  an  infinite  self-velocity.  The  discrete  system 
is  developed  into  a  lagrangian  numerical  method  which  exactly  preserves  the  kinetic  energy, 
impulse  and  angular  momentum  of  the  incompressible  fluid  flow.  We  compare  the  method  with 
standard  vortex  methods. 

We  also  present  results  of  a  finite  difference  method  based  on  the  magnetization  formulation  of 
incompressible  fluid  flows.  The  magnetization  formulation  has  advantages  over  a  formulation  in 
primitive  variables  in  that  the  pressure  does  not  appear  explicitly  in  the  equation  of  motion  for  the 
magnetization  variable.  We  extend  the  continuous  equations  for  the  case  of  variable  density  flows 
and  discuss  the  extension  of  the  methods  to  the  variable  density  case. 


TURBULENT  EDDY  STRUCTURES, 
COMBUSTION  AND  CHEMICAL  REACTIONS 


Jean-Pierre  CHOLLET 
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Turbulence  strongly  affects  transport  and  dispersion  of  contaminants 
advected  by  velocity  fields.  For  a  fluid  made  of  species  which  undergo 
chemical  reactions,  concentration  fields  evolve  with  time  and  space  under 
competing  effects  of  flow  turbulence,  chemical  reaction  and  diffusion 
fluxes . 

The  mixing  layer  has  been  selected  as  a  prototype  of  shear  flows  which 
develop  vortex  structures  from  unstability  mechanisms.  Such  vortices  are 
observed  to  strongly  affect  mixing  and  then  chemical  activity.  First, 
binary  reactions  are  considered  in  numerical  simulations  of  three 
dimensional  incompressible  mixing  layers  using  pseudo  spectral  methods  with 
periodic  boundary  conditions  (i.e.  temporal  approximation).  Reaction  rate 
is  observed  to  be  the  quantity  of  interest  to  study  chemical  activity  and 
to  trace  some  properties  of  turbulent  fields,  especially  in  the  strained 
regions.  Chemical  activity  is  observed  to  develop  along  vortex  structures, 
with  increasing  spatial  intermittency  as  the  reaction  gets  faster  (higher 
Damkholer  number).  Depending  on  the  way  vorticity  fields  are  triggered 
initially,  resulting  species  concentrations  evolve  differentely  with  more 
or  less  three  dimensionality.  Results  from  a  preliminary  lagrangian 
particle  tracing  agree  with  this  general  trend.  Chemical  reactions 
emphasize  also  the  contribution  from  diffusivity  and  viscosity  (Schmidt 
numbers)  and  then  enlights  a  discussion  about  direct  versus  large  eddy 
simulations  and  the  need  for  subgrid  scale  models. 

Combustion  differs  from  other  reacting  flows  because  of  density 
variations  induced  by  strong  heat  release.  Chemistry  is  complex  and 
temperature  dependent.  A  numerical  method  has  been  derived  specifically  in 
order  to  compute  evolutions  in  time  and  space,  even  with  shocks  and  steep 
temperature  or  concentration  gradients.  Wide  ranges  of  flow  velocities  have 
been  considered  up  to  supersonic  conditions.  Flow  structure  gets  rather 
intricate  with  both  evolving  vortex  structures  and  non  stationary  shock 
patterns.  Simulations  have  been  developed  in  two  problems  of  interest  for 
aircraft  engine  design,  the  underexpanded  jet  and  the  mixing  layer  confined 
between  rigid  walls.  Results  from  simulations  agree  with  those  classicaly 
derived  from  instability  theories.  Moreover  simulation  provides  pictures  of 
time  and  space  distribution  of  vorticity  along  the  flow.  Multiple  step 
reactions  have  been  also  considered  in  the  simulations,  assuming  either 
species  partial  equilibrium  or  explicit  evolution  of  species  mass 
fractions,  depending  on  the  ratios  of  charatestic  chemical  times  to 
turbulence  time  scale. 
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VORTEX  GENERATION  AND  EVOLUTION  IN 
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The  experiments  of  mixing  layers  performed  at  CalTech  by  Roshko  and  his  group  [1]  showed 
the  presence  of  large-scale  organized  vortices  even  extremely  far  downstream  of  their  origin.  It 
turned  out  that  these  vortices  resulted  from  the  growth  of  hydrodynamic  instabilities  yielding 
roll-up  of  fundamental  vortices  and  successive  pairings  of  the  latter.  These  trends  were  retrieved 
numerically  in  particular  by  Normand  et  al.  [2],  who  developed  two-dimensional  direct  numerical 
simulations  of  a  spatially-growing  mixing  layer  forced  upstream  by  a  small  random  perturbation. 
However  all  natural  plane  shear  layers  are  known  to  develop  three-dimensionality,  both  in  the 
large  and  small  scales,  characterized  in  particular  by  thin  streamwise  hairpin  vortices  streched 
between  the  large-scale  rollers  [3]  and  developed  turbulence  with  k~ S/<3  Kolmogorov  spectra. 
We  will  present  three-dimensional  simulations  of  both  incompressible  and  compressible  mixing 
layers  reproducing  these  features.  They  also  point  out  a  great  sensitivity  of  the  flow  structure  to 
initial  conditions:  in  some  cases,  the  flow  pattern  appears  to  be  controlled  by  unstable  oblique 
modes  which  are  not  the  most  amplified  ones  predicted  by  the  linear  theory  of  stability  [4]. 

Nevertheless  this  does  not  necessarily  contradicts  a  possible  universality  of  transitional 
mechanisms  leading  to  developed  turbulence:  we  will  also  show  results  of  simulations  of  jets 
and  boundary  layers  (both  in  incompressible  and  compressible  situations)  in  which  analogous 
instability  mechanisms  are  at  work,  yielding  similar  vortical  structures  (at  least  to  the  lowest 
order  of  approximation);  in  particular,  the  staggered  or  aligned  A-vortices  found  in  boundary 
layers  close  to  the  wall  also  result  from  the  growth  of  oblique  modes  [5]. 

We  will  also  present  comparisons  of  results  obtained  by  means  of  different  numerical  meth¬ 
ods,  in  particular  the  vortex  method  used  by  Ashurst  and  Meiburg  [6]  which  were  found  to 
reproduce  satisfactorily  the  coherent  structures  of  incompressible  mixing  layers,  free  jets  and 
wakes  during  the  early  stages  of  transition.  Nevertheless,  when  turbulence  develops  in  the  small 
scales,  dissipative  mechanisms  have  to  be  taken  into  account:  the  molecular  diffusion,  which  will 
act  at  scales  smaller  than  the  mesh  if  the  Reynolds  number  is  high  enough.  In  this  case,  it  has 
to  be  modelled  by  means  of  appropriate  techniques  (subgrid-scale  parameterization),  which  will 
be  briefly  discussed. 
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VORTICITY  ON  THE  BOUNDARIES  FOR  THE  2D  NAVIER  STOKES 
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Some  results  concerning  the  equivalence  between  the  Navier-Stokes  equations  in  primitive 
variables  on  one  hand  and  in  velocity  voticity  formulation  on  the  oher  hand,  are  recalled. 

To  compute  the  velocities,  we  made  the  choice  to  use  Poisson  equations.  Thus,  we 
present  the  conditions  that  the  vorticity  u>  must  satisfy  on  solid  boundaries  in  order  to 
enforce  the  definition  of  u  as  the  curl  of  the  velocity,  and  therefore  the  mass  conservation. 

In  the  unsteady  computations  that  are  presented, (driven  cavity,  flow  around  a  circular 
cylinder),  these  conditions  are  satisfied  by  means  of  an  influence  matrix  technique  that 
ensures  at  each  time  step ,  the  definition  of  u  and  the  vanishing  of  the  divergence  of  the 
velocity  field. 


Simulation  of  2D  Turbulence  at  near-infinite  Reynolds  number 


David  G.  Dritschel 


The  “moment  accelerated  contour  surgery”  (MACS)  algorithm  has  made  possible  an 
altogether  new  approach  to  2D  turbulence.  It  has  enabled  for  the  first  time  the  simulation 
of  2D  vortex  flows  at  Reynolds  numbers  enormously  higher  than  was  previously  possible. 

MACS  is  an  extended  “contour  dynamics”  method  based  on  the  Lagrangian  descrip¬ 
tion  of  fluid  mechanics.  Material  conservation  of  vorticity,  the  fundamental  dynamical 
property  of  inviscid  fluid  motion,  is  satisfied  automatically.  MACS  builds  on  contour 
dynamics  by  combining  it  with  other  techniques  in  a  new  way  that  results  in  a  ten-  to 
hundred-fold  gain  in  computational  efficiency. 

The  results  of  three  large  ensembles  of  calculations  differing  only  in  spatial  resolution 
demonstrate  rapid  convergence  to  essentially  inviscid  dynamical  behavior.  The  statistical 
vortex  properties  obtained  by  and  large  differ  substantially  from  all  previously-reported 
results  for  two-dimensional  turbulence.  In  particular,  many  more  small-scale  vortices  are 
observed,  with  the  number  density  distribution  not  following  a  power  law  but  rising  pro¬ 
gressively  more  steeply  with  decreasing  vortex  size.  Also,  and  as  a  consequence,  the  time 
decay  of  basic  flow  properties  (e.g.,  total  enstrophy)  is  significantly  slower,  and,  moreover, 
does  not  fit  the  recently-proposed  universal  scaling  theory  for  two-dimensional  turbulence. 

The  ability  of  MACS  to  model  a  greatly  extended  range  of  spatial  scales,  as  compared 
to  conventional  models,  and  to  preserve  vorticity  gradients  from  artificial  numerical  erosion, 
is  what  distinguishes  its  performance  and  potential  from  that  of  conventional  Eulerian 
methods. 


Computer  simulation  of  fluid  flow  with  heat  and  mass  transfer  in 

engineering  and  science 
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Industrial  applications  require  simple  and  unified  computer  simulation  of  different  rather  complex 
problems  performed  by  users  from  common  research  and  development  staff  without  special 
training  in  numerical  mathematics.  Two  multipurpose  Applied  Code  Packages  of  Computational 
Fluid  Dynamics  named  AEOL  and  NEPTUNE  involve  multi-dimensional  models  (the  Euler 
equations  for  compressible  gas  and  the  Navier-Stokes  equations  for  incompressible  fluid)  with 
extensive  physics  and  chemistry.  They  may  be  easily  appended  and  modified  by  a  user  in  a 
straightforward  manner  and  may  be  incorporated  in  Computer  Aided  Design  systems  or  used  for 
student  training.  The  main  distinction,  which  made  them  particularly  suitable  in  engineering,  is 
inplementation  of  the  conservative  homogeneous  schemes  on  rectangular  grids  in  the  physical 
domain,  allowing  flow  simulation  in  regions  of  complicated  shape,  even  varying  in  time.  The 
simplicity  of  algorithms  provides  implementation  on  transputor  and  multiprocessor  systems,  vector 
and  array  processors,  etc  .. 


The  code  performances  on  PC  may  be  demonstrated. 


N  W^pects  of  Vortex  Dynamics  by  Numerical  Simulation: 

_ cal  Wave  Dynamics,  and  Vortex-Fine  Scale  Turbulenc* 

by  Fazle  Hussain 

By  investigating  a  straight  axisymmetric  tube  having  an  initially  nonunifonn  core  size 
along  its  axis,  with  and  without  superimposed  fine-grained  turbulence,  we  bring  new  life  to 
local  heiicity  concepts.  While  the  relative  helicity  and  heliciry  densities  are  of  little  use,  we  find 
that  the  separation  of  the  flow  field  into  right  and  left  handed  components  by  the  "complex  helical 
wave  decomposition"  is  highly  useful.  This  decomposition  seems  to  open  up  an  entirely  new 
discipline  of  truly  3D  vortex  dynamics  (including  core  and  wave  packet  dynamics),  as  well  as  to 
provide  a  rich  new  set  of  tools  for  analyzing  3D  vortical  fields  with  coherent  structures.  Wc  find 
that  the  evolution  of  the  laminar  vortex  can  be  explained  both  by  traditional  means  (such  as 
coupling  between  swirling  and  meridional  flow)  and  in  terms  of  helical  structures;  in  the  context 
of  such  structures  the  dynamics  is  much  simpler,  more  revealing,  and  allows  better  qualitative 
predictions  of  the  flow’s  evolution  than  docs  the  well-known  vortex  dynamics.  We  show  that 
vortex  core  dynamics,  ignored  by  all  texts  and  most  researchers,  can  have  a  significant  effect  on 
a  vortex's  evolution  and  interaction  with  others  or  with  fine-scale  turbulence.  Core  dynamics  is 
better  understood  in  terms  of  nonlinearly  interacting  polarized  wave  packets  whose  evolution 
equation  is  extracted  from  the  vorticity  transport  equation  by  use  of  projection  operators.  The 
core  size  oscillation  is  controlled  by  viscous  effects  (decaying  slower  at  high  Re),  the  oscillation 
frequency  decreasing  with  increasing  Re,  but  having  a  finite  limit  as  Re-*«\  Internal  dynamics 
disappears  first  near  the  outer  edge,  gradually  retracting  to  an  increasingly  smaller  region  around 
the  axis. 

For  the  turbulent  vortex  we  find  the  striking  phenomenon  of  organization  of  incoherent 
vorticity  in  the  boundary  layer  surrounding  the  vortex  into  spiral  type  structures  that  produce 
intennittency,  and  the  inviscid  generation  of  increasingly  larger  scales  -  a  prime  example  of  non¬ 
trivial  coupling  between  large  and  small  scales  in  turbulent  flows.  Our  finding  is  in  contrast  with 
the  prior  suggestion  of  formation  of  sheet- like  intermittency  structures,  and  we  find  the 
azimuthally  polarized  smaller-scale  vortices  to  anti-cascade  into  larger  structures  by  pairing.  The 
organization  is  accompanied  by  a  strong  tendency  for  the  right  and  left  handed  components  of  the 
vorticity  field  to  separate  spatially.  Moreover,  the  small  scales  take  energy  form  the  coherent 
vortex  and  can,  if  the  Reynolds  number  is  sufficiently  high,  feedback  on  the  coherent  vortex  in 
such  a  way  as  to  excite  non-axisymmetric  bending  waves.  This  type  of  feedback  or  backscatter 
as  a  result  of  coupling  between  large  and  small  scales  can  not  be  molded  by  any  type  of  eddy 
viscosity.  Our  study  raises  some  doubts  about  the  validity  of  the  local  isotropy  concept  to 
turbulent  shear  flows.  We  conclude  with  a  cascade  model  for  turbulence  scale  hierarchy. 


Vortex  Method  Approxmations  for  Two-Phase  Flows  with  Surface  Tension 


Thomas  Y.  Hou 
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A  variety  of  physical  phenomena  involve  propagating  interfaces.  It  ranges  from  the 
familiar  water  wave  problem  to  problems  in  flame  propagation  and  dendritic  solidification. 
It  is  well-known  that  the  vortex  sheet  problem  is  ill-posed  and  singularities  are  formed  in  a 
finite  time  without  physical  regularizations.  Here  we  consider  the  stabilizing  effect  of  surface 
tension  for  interfacial  flows,  including  water  waves,  Hele-Shaws,  and  two-phase  flows.  The 
stabilizing  effect  of  surface  tension  was  not  well  understood  before,  partially  due  to  the 
difficulty  in  distiguishing  the  numerical  instability  from  a  phyiscal  one.  In  our  study,  we 
use  a  spectrally  accurate  point  vortex  method  to  approximate  the  governing  equations.  We 
show  that  this  method  is  stable  as  long  as  the  physical  solution  is  smooth.  In  the  process  of 
proving  this  result,  we  develop  a  powerful  analytical  technique  to  analyze  stability  of  vortex 
methods  by  a  subtle  energy  estimate  without  using  analyticity.  Our  study  also  indicates 
that  there  is  a  strong  compatibility  constraint  that  must  be  satisfied  by  a  numerical  method 
in  order  to  be  stable.  Violation  of  this  constraint  will  lead  to  numerical  instabilities.  For 
example,  if  one  uses  an  alternating  trapezoidal  rule  to  approximate  the  velocity  integral 
and  uses  a  cubic  spline  approximation  to  the  curvature,  the  resulting  scheme  is  numerically 
unstable. 

Another  aspect  of  our  study  is  to  investigate  the  question  of  whether  or  not  surface 
tension  can  prevent  singularity  formation.  This  has  been  an  unsettled  question  for  a  long 
time.  Since  our  method  is  proved  to  be  stable  and  spectrally  accurate  as  long  as  the  solution 
is  smooth,  this  provides  a  very  accurate  method  to  compute  the  singularity  if  there  is  one. 
Our  study  shows  that  if  surface  tension  is  above  some  critical  value,  the  interface  problem 
has  a  smooth  solution  for  all  time.  If  surface  tension  is  below  the  critical  value,  the  curvature 
of  the  interface  grows  exponentially  in  time.  This  makes  it  extremely  difficult  to  distinguish 
numerically  the  exponential  growth  from  a  finite  time  singularity. 

One  major  difficulty  in  computations  of  interfacial  flows  with  surface  tension  is  the  se¬ 
vere  time  step  constraint;  e.g.  At  must  be  less  than  h3  for  an  explicit  time  discretization.  In 
our  study,  we  propose  to  use  an  implicit  time  discretization  coupled  with  an  efficient  precon¬ 
ditioned  iterative  method.  This  combination  provides  a  very  efficient  numerical  algorithm 
-  a  factor  of  0(JV2)  faster  than  the  usual  explicit  discretization.  Here  the  choice  of  the 
preconditioner  is  crucial.  This  is  constructed  from  the  leading  order  solution  operator  in  its 
normal  and  tangential  coordinates,  guided  from  our  stability  analysis.  Extensive  numerical 
experiments  will  be  presented. 


"A  COUPLED  POTENTIAL-BOUNDARY  LAYER  CALCULATION 
METHOD  FOR  UNSTEADY  FLOWS  AROUND  AIRFOILS" 

M.  Kermarec,  A.F.  Decaix,  P.  Renon,  D.  Favier,  C.  Maresca 
IM2-IMFM-UM34  CNRS,  163  Av.  de  Luminy,  13009  Marseille 

A  numerical  method,  based  on  coupling  a  potential  model,  and  a 
boundary-layer  calculation,  has  been  developed  to  predict  unsteady  airloads  of 
airfoils  oscillating  in  two-dimensional  flows.  One  original  feature  of  the 
numerical  approach  consists  in  determining  the  geometry  of  the  boundary  at  the 
junction  between  the  two  calculation  domains.  The  coupling  boundary  location  is 
specifically  determined  for  flow  zones  including  boundary-layer  separation  on 
the  airfoil  surface  and  downstream  of  the  trailing-edge,  in  order  to  verify  the 
following  conditions  :  the  coupling  boundary  represents  a  streamline  in  the 
frame  linked  to  the  oscillating  airfoil;  a  constant  pressure  value  is  imposed  on 
the  separated  boundary  layer  coupling  line,  equal  to  the  value  at  the  trailing 
edge;  pressure  values  across  the  wake  behind  the  trailing  edge  must  be  equal. 

The  numerical  approach  has  first  been  validated  in  steady  flow 
configurations  around  the  model  at  rest.  Comparisons  between  calculation  and 
experiment  have  been  performed  on  overall  lift  and  drag  coefficients.  A  new 
transition  criterion  has  been  deduced  from  boundary-layer  velocity  profiles 
measured  on  a  flat  plate  model.  Using  this  transition  criterion  within  the 
numerical  approach,  a  good  agreement  between  calculations  and  experiments  has 
been  obtained  on  lift  and  drag  coefficients  of  the  airfoil  including  different 
airfoil  geometry  (  NACA0012  and  ONER  A  OA209  )  and  Reynolds  numbers 
varying  from  10^  to  3x10^. 

For  unsteady  flow  configurations,  the  first  step  of  the  numerical 
approach  has  been  to  introduce  a  modification  of  the  classical  mixing  length 
turbulence  model.  This  modified  turbulence  model  has  been  based  on  physical 
arguments,  and  also  deduced  from  velocity  measurements  performed  within  the 
boundary-layer  of  the  flat  plate  model  oscillating  either  in  fore  and  aft  motion 
or  in  pitching  motion.  The  numerical  method  has  been  applied  to  the  pitching 
airfoil  case  below  dynamic  stall,  at  low  and  moderated  amplitudes  of  incidence 
oscillations.  The  calculated  airloads  coefficients  are  shown  to  well  match  the 
experimental  hysteresis  loops.  Moreover,  the  comparisons  between  calculations 
and  experiments  also  exhibited  the  strong  dependence  of  the  numerical  results  on 
the  good  prediction  of  the  geometrical  boundary  location  between  the  two 
calculation  domains. 

Finally,  the  present  numerical  coupling  method  has  been  checked  for 
airfoils  oscillating  in  pitch  through  dynamic  stall.  In  this  case,  the  boundary- 
layer  calculation  is  not  performed,  and  the  unsteady  motion  of  the  turbulent 
separation  point  is  deduced  from  local  surface  measurements  on  the  upper  side 
of  the  airfoil.  Potential  flow  results  obtained  with  such  experimental  data,  and 
with  the  numerical  determination  procedure  of  the  coupling  boundary  location, 
show  the  capability  of  the  method  to  predict  the  experimental  lift  coefficient 
behavior  through  stall. 
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ABSTRACT 

The  linear  instability  of  a  family  of  inviscid,  two-dimensional,  variable-density  shear  layers 
and  wakes  is  investigated.  Vorticity  profiles  corresponding  to  a  monotonically  increasing  velocity 
profile  are  first  examined.  A  larger  family  of  initial  vorticity  distributions  which  model  the  merger 
of  two  unequal  vorticity  layers  of  opposite  sign  is  then  considered.  The  latter  is  obtained  by 
superimposing  on  the  former  a  wake  component,  characterized  by  a  spread,  8,  and  a  velocity 
deficit,  W.  The  initial  density  distribution  resembles  a  temperature  spike  and- is  described  by  a 
lui^kness,  o,  and  a  temperature  ratio,  Tr.  The  stability  properties  of  the  layers  are  interpreted  in 
terms  of  a  four-dimensional  parameter  space  (W,8,Tr,o).  Finally,  the  non-linear  evolution  of  the 
flowfieid  is  illustrated  using  the  transport  element  method. 

Flowficld  stability  exhibits  strong  sensitivity  tc  the  details  of  the  density  distribution.  In 
the  absence  of  the  wake  ->mponent,  the  stability  properties  of  the  heated  layer  are  divided  into 
three  categories  according  to  the  thickness  of  the  density  profile,  a,  and  the  vorticity  thickness,  8w. 
For  <y  »  8w,  instability  of  •:  Kelvin-Helmholtz  mod?  in  a  uniform-density  flow  is  recovered. 
When  o  —  8\v,  the  shear  mode  is  inhibited,  whi,e  this  trend  persists  for  o  ^  6^,  the  layer 
becomes  characterized  by  the  appearance  o i  additional  short-wavelength  unstable  modes  which 
become  dominant  as  o  decreases  and  Tr  increases.  Addition  of  a  wake  component  is  shown  to  alter 
this  behavior,  and  to  oppose  the  stabilizing  effects  of  heat  release.  In  this  case,  the  shear  layer 
mode  always  dominates  the  wake  mode,  and  the  presence  of  heated  sublayer  has  a  weak  effect  on 
the  instability  of  the  layer  when  8  is  large,  but  may  influence  the  phase  speed  of  unstable  waves 
whenever  the  zones  of  high  vorticity  and  high  density  gradient  coincide. 


Large  Scale  Simulations  Using  Vortex  Methods. 
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Abstract. 

Simulations  are  presented  for  bounded  and  unbounded  viscous  incompressible 
flows.  Our  numerical  method  is  based  on  vortex  methods.  The  classical  invh.'ud 
scheme  is  enhanced  to  account  for  viscous  effects  via  the  method  of  particle  strength 
exchange.  The  method  is  extended  to  account  for  the  enforcement  of  the  no-slip 
boundary  condition  as  well  by  appropriately  modifying  the  strength  of  the  parti¬ 
cles.  Computations  are  possible  for  extended  times  by  periodically  remeshing  the 
vorticity  field. 

XL*  particles  are  advanced  using  the  Biot-Savart  law  for  the  evaluation  of  the 
velocity.  Computations  are  made  using  C?(105)  vortex  particles  by  efficiently  imple¬ 
menting  the  method  of  multipole  expansions  for  v*  •  ?  .-  computer  architectures. 

The  method  is  used  to  simulate  the  merging  of  Gaussian  vortic. 
bounded  fluid  as  well  as  unsteady  flows  behind  circular  cylinders  for  a  v. .  range 
of  Reynolds  numbers.  Direct  comparisons  are  made  with  spectral  method  for  the 
unbounded  flows  and  with  a  variety  of  numerical  methods  and  experimental  result  - 
for  the  bounded  flows. 
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Vortex  Sheet  Computations 

Robert  Krasny 
University  of  Michigan 

I  will  discuss  recent  work  related  to  the  computation  of  vortex  sheet 
motion.  The  following  problems  will  be  addressed: 

1.  Periodic  vortex  sheet  roll-up,  convergence  of  the  vortex  blob  method  past 
the  critical  time,  structure  of  the  rolled-up  spiral. 

2.  Wake  patterns  computed  by  the  vortex  blob  method  and  comparison  with 
experimental  flow  visualization  by  Couder  and  Basdevant  (JFM,  vol.  173, 
1986). 

3.  Vortex  sheet  separation  at  a  sharp  edge,  comparison  with  Pullin’s  com¬ 
putation  of  self-similar  vortex  sheet  roll  up  past  a  semi-infinite  flat  plate 
(JFM,  vol.  88,  1978),  jet  formation  and  instability,  vortex  pair  formation 
due  impinging  flat  plates. 


VORTEX  BLOB  SIMULATION  OF  TWO-DIMENSIONAL  FLOW 


IN  CHANNEL  WITH  COMPLEX  GEOMETRIES. 

Henryk  Kudela 

Technical  University  of  Wroclaw 
50-370  Wroclaw,  Wybrzeze  Wyspiahskiego  27,  Poland 

Mostly  the  vortex  blob  method,  originally  proposed  by  A.J.  Chorin,  was  used  to  numerical  solution 
of  the  Navier-Stokes  equations  at  high  (as  well  as  at  low  )  Reynolds  number  in  area  with  relatively 
simply  geometry  (a  backward-facing  step,  a  cylinder,  a  cavity).  To  fulfill  the  condition  of  cancellation 
of  the  normal  velocity  to  the  wall  one  has  to  add  to  the  velocity  generated  by  the  blobs  the  auxiliary 
potential  flow.  For  such  simply  area  as  mentioned  above,  that  auxiliary  potential  flow  can  easily 
calculated.  In  order  to  be  able  to  apply  the  vortex  blob  method  to  any  geometry  (in  two  dimension) 
we  used  the  fast  elliptic  solvers  together  with  the  capacitance  matrix  technique.  The  solution  of  that 
boundary  value  problem  for  this  auxiliary  potential  flow  does  not  depend  on  any  gradients  of  the 
flow.  Potential  velocity  of  each  blobs  were  calculated  by  area-weighting  scheme.  The  tangential 
boundary  condition  is  realized  through  the  generation  of  the  vortex  blob  at  the  wall.  The  generations 
was  carried  out  at  points  spaced  at  a  distance  h  along  the  wall.  We  used  h=0.1,  and  a=h^-^5  where 
a  is  cut-off  radius.  Ordinary  differential  equations  for  displacements  of  the  blobs  were  solved  by  the 
first  order  Euler  methods  with  time  step  At=0.05.  Numerical  grid  for  elliptic  solvers  was  100*20. 
The  ratio  of  the  length  to  height  of  the  channel  was  as  one  to  ten.  We  did  not  use  the  vortex  sheet 
method  for  flow  near  the  boundary.  To  test  and  verify  of  the  program  we  studied  the  flow  over 
backward-facing  step  for  different  Reynolds  numbers  (see  Ghoniem,Cagnon  (1987),  and  Sethian, 
Ghoniem  (1988),  J.Comp.Phys.)  and  we  obtained  good  agreement  with  those  work.  We  carried  out 
also  calculation  for  flow  over  rectangular  step  and  for  flat  plate  normal  to  the  flow  in  channel.  The 
results  (velocity  profiles,  streamlines,  instantaneous  portrait  of  the  vortex  blobs)  were  promising. 


Particle  Approximation  of  the  Vlasov-Poisson  Equation 
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In  this  work  we  consider  the  spatially  periodic  1-D  single  component  Vlasov-Poisson  equation 
from  plasma  physics.  The  initial  condition  is  the  analogue  of  a  vortex  sheet  which  we  call  an 
electron  sheet.  It  consists  of  a  delta  function  supported  over  a  curve  in  x-v  space.  A  number  of 
exact  solutions  to  this  problem  are  known.  Some  of  them  exist  for  all  time,  while  others  have  a 
finite  time  after  which  the  solution  of  the  Vlasov-Poisson  equation  develops  singularities. 

We  approximate  the  Vlasov-Poisson  equation  with  electron  sheet  initial  data  by  a  particle  method. 
We  demonstrate  that  the  solution  of  the  particle  method  converges  to  the  exact  solution  over  its  time 
interval  of  existence.  We  discuss  qualitative  properties  of  the  solution  after  the  formation  of  a 
singularity. 


Point  Vortices  and  Localization  in  Euler  Flows. 
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It  is  very  old  and  natural  the  attempt  to  study  an  incompressible  non-viscous  fluid 
when  the  vorticity  in  concentrated  in  a  subset  of  physical  space  of  low  dimensionality  .  For 
planar  situation  die  first  study  of  a  fluid  when  the  vorticity  is  sharply  concentrated  in  very 
small  regions  (schematized  as* points)  has  been  done  in  the  last  century,  introducing  the  point 
vortex  system .  A  formal  connection  between  this  model  and  the  Euler  equation  is  discussed 
in  many  textbooks.  However  only  recently  a  rigorous  connection  have  been  proved.  A  short 
r  eview  of  this  topics  is  the  argument  of  the  present  talk.  In  particular  it  is  possible  to  prove  a 
>c.  n  of  "localization"  of  the  motion  in  the  sense  that  the  evolved  vorticity  is  concentrated  in  N 

regions  of  diameter  d  which  vanishes  as  e  (the  initial  diameter )  vanishes  for  any  fixed  time. 
For  smooth  problems  this  is  a  direct  consequence  of  the  continuity  with  respect  to  the  initia 
datum.  But  our  problem  is  singular  because  the  initial  datum  is  sharply  concentrated  and  the 

time  evolution  (via  the  Euler  Equation)  has  a  logarithmic  divergent  kernel.  Actually,  when  e  is 
very  small,  the  field  in  each  blob  becomes  very  large  and  it  is  difficult  to  exclude  that  thin 
filaments  of  vorticity  are  pushed  away.  We  can  prove  that  this  fact  does  not  happen.  As  main 
consequence  we  can  prove  a  rigorous  connection  between  the  Euler  equation  and  the  point 
vortex  theory. 

Some  results  in  this  direction  has  been  obtained  in  the  last  ten  years  but  only  recently 
has  been  completely  solved  [1]. 

More  exactly,  consider  an  initial  datura  of  the  form: 

N 

<Oe(i>0)  =  <Oe;i(2L»0)  (1) 

where  is  a  function  with  a  definite  sign  supported  in  a  region  Ae;j  CR2  suCh  that 

Ae;i  s  supp  <s%;i  C  !  e)  ;  I(zj  I  e)  fi  I  e)  if  i*j  (2) 

for  e  small  enough  and  where  X(i,r)  is  a  circle  of  center  z  and  radius  r . 

Moreover 

|  di  We;j(i,0)  =  ai  e  3R  ( vortex  intensity )  (3) 


and 


i  (Oe;i(&,0) !  5  const  £*h 


(4) 


We  can  prove  the  following  result: 

Blfiggffl 

Denote  by  o>e(jc,t)  the  time  evolution  of  oicfoO)  via  the  Euler  Equation.  Then 


i)  for  all  d  >  0  there  exists  Co(d,T)  such  that,  if  £  <  £o(d,T) ,  then 

supp  C0e;i(2,t)  C  I(2i(t)  i  d)  for  any  te  [0,T]  (5) 

where  gi(t)  is  the  solution  of  the  ordinary  equation  ijaieni  ( called  point  vertex  system) 


3ra(t>  -  -  i  ii(t)  -  Dio 


2i(0)  =  ii 


(6) 


ii)  for  any  continuous  bounded  function  f(a) 


^  J  di  coc(2i.,t)  f(i)  =  1  ai  *Gu(0) 

e — >0  i=l 


(7) 


We  shortly  comment  the  statement  of  the  theorem.  Position  i)  states  that  the  blobs  of 

vorticity  remain  localized  until  time  T  for  any  d  and  T  provide  we  choose  e  small  enough. 
Position  ii)  states  that 

weak  jsf 

cae(s,t) - *  X  ai  *<*1(0)  (8) 

e-»0  t=l 

where  5()  denotes  the  Dirac  measure.  This  last  statement  gives  a  rigorous  justification  of  the 
point  vortex  model. 

In  three  dimension  the  situation  is  less  satisfactory.  The  naive  generalization  of  point 
vortices  would  be  filaments  of  vorticity.  However  a  simple  calculation  show  that  they  move 
(in  general)  with  an  infinite  speed  and  so  a  Theorem  like  the  previous  one  is  hopeless.  Some 
point  systems,  called  vortons,  has  been  studied,  but  they  do  not  conserve  some  important 
quantities  and  so  are  less  fundamental  than  the  point  vortex  in  two  dimension.  However  the 
can  be  important  for  numerical  purpose. 

For  a  review  of  the  topics  discussed  in  this  talk,  see  [2], 

1.  Marchioro  C.,  Pulviremi  M.,:  Vortices  and  localization  in  Euler  flows,  Commun,  Math. 
Phys.  ( sumbitted  to) 

2.  Marchioro  C.,  Pulvirenti  M.,  Mathematical  Theory'  of  Incompressible  Non-Viscous 
Fluids,  Berlin,  Heideiberg, New  York  Springer  Verlag  (in  press)  (1992) 


A  PARTICLE  IN  CELL  METHOD  FOR  THE  COMPRESSIBLE  EULER  EQUATIONS 

S.  Mas-Gallic  (*),  M.  Louaked  (**) 


SUMMARY 

The  flow  of  an  inviscid  compressible  fluid  is  governed  by  the  compressible  Euler  equations, 

757  +  (u.V)u  +  -  Vp  =  0, 

a t  p 

^  +  V(pu)  =  0, 

P  =  fcp7, 

where  uf  p  and  p  are  respectively  the  velocity,  the  density  and  the  ptessure  of  the  fluid  and  where  7  is  the 
specific  heat  ratio.  We  consider  the  two  dimensional  case  and  introduce  the  Helmholtz  decomposition  of  the 
velocity  field  u 

u  =  V<t>  -I-  v 

where  v  is  the  rotationnal  part  of  the  velocity.  Introducing  the  vorticity  w  =  Vxu=Vxv,  the  system  of 
Euler  equations  is  then  written  under  the  following  form 

^  +  V(p(V«i  +  v))  =  0, 

^  +  V(W(V*  +  v))  =  0, 

-Av  =  V  x  w 

^  +  i|  v  +  v^  +  tviI^-  =  ,, 

A q  =  w2  +  Vu>  x  (V<j>  +  v). 

Under  this  form  we  notice  that  both  the  density  p  and  the  vorticity  u>  are  convected  by  the  flow  with  the 
velocity  u  =  V^+ v.  Then  we  consider  the  equation  which  gives  the  evolution  of  the  potential  <j>.  Introducing 
the  hamiltonian 


this  equation  is  rewritten 


H(xP)  =  -W\ 


^  +  v.Vd  =  q  -  ^|v|2  - 


and  we  recognize  a  Hamilton- Jacobi  equation  with  a  convection  term. 

The  method  of  resolution  of  the  Euler  compressible  system  is  based  on  these  two  remarks.  The  convection 
equations  are  solved  by  a  particle  method  whereas  the  Hamilton-Jacobi  equation  is  solved  by  the  finite 
difference  scheme  proposed  by  S.  Osher  and  J.  Sethian  and  the  Laplace  equation  by  a  classical  finite  difference 
scheme.  A  grid-particle  operatorand  a  particle-grid  operator  are  then  defined,  as  in  any  particle  in  cell 
method.  The  grid-particle  operator  associates  a  particle  approximation  to  the  quantities  defined  on  the  grid, 
the  particle  velocities  are  computed  with  this  procedure,  whereas  the  particle-grid  operator  defines  a  grid 
approximation  for  the  quantities  associated  to  the  particles.  In  order  to  ensure  the  conservativity  of  the 
scheme,  the  particle-grid  operator  is  defined  following  the  ideas  of  J.U.  Brackbill  and  H.M.  Ruppel. 


(*)  Centre  de  Mathematiques  Appliquees,  Ecole  Poly  technique,  91128  Palaiseau  Cedex,  France 

(**)  Laboratoire  d’Analyse  Numerique,  T55-65,  5e  etage,  Universite  Pierre  et  Marie  Curie,  4  Place  Jussieu, 
75252  Paris  Cedex  05,  France 
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We  study  the  inviscid  mechanisms  governing  the  three  dimensional  evolution  of  round  jets  by 
means  of  vortex  dynamics  simulations.  A  typical  simulation  starts  from  an  unperturbed 
axisymmetric  layer  of  vortex  filaments.  For  non-swirling  jets,  the  filaments  initially  have  the  form 
of  rings,  while  helical  vortex  filaments  introduce  swirl  into  the  jet.  The  axisymmetric  jet  shear  layer 
can  then  be  perturbed  in  different  ways  by  modulating  the  strength  of  the  filaments  (corresponding 
to  acoustic  forcing  in  the  axial  direction)  or  by  displacing  the  filament  centerlines  (thus  simulating 
the  effect  of  corrugated  or  indented  nozzle  shapes  or  acoustic  forcing  in  the  azimuthal  direction). 
The  spatially  periodic  calculations  provide  a  detailed  picture  of  the  processes  leading  to  the 
concentration,  reorientation,  and  stretching  of  the  vorticity  for  different  combinations  of 
perturbations.  In  the  purely  axisymmetric  case,  a  wavy  perturbation  in  the  streamwise  direction 
leads  to  the  formation  of  vortex  rings  connected  with  braid  regions,  which  become  depleted  of 
vorticity.  An  additional  azimuthal  wave  leads  to  a  fully  three-dimensional  evolution,  which  results 
in  concentrated  braid  vortices  as  well  as  a  potential  vortex  ring  instability.  For  non-swirling  jets 
subject  to  axisymmetric  and  azimuthal  waves,  we  observe  the  formation  of  counterrotating 
streamwise  braid  vortex  pairs,  while  the  presence  of  helical  waves  or  swirl  in  the  jet  leads  to  the 
formation  of  braid  vortices  all  of  the  same  sign. 
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Abstract 


Vortex  generation  in  a  rotating  tank  by  a  convergence  of  vortidty  caused 
by  suction  through  a  tube  centered  on  the  rotation  axis  is  reconsidered  in  view  of 
an  analysis  of  Ekman  layers  dynamics.  New  experimental  studies,  based  on  flow 
visualisations  and  velocity  profile  measurements  by  Laser  Doppler  Anemometry, 
are  presented,  which  show  that  the  whole  secondary  flow  is  passing  through  the 
Ekman  layers.  Ekman  layers  therefore  govern  the  structure  of  the  flow.  A  model, 
inspired  by  Stem  (1975),  establishes  that  the  circulation  of  the  vortex  is  equal  to  the 
flowrate  divided  by  the  Ekman  layer  depth.  When  the  Ekman- number,  build  on 
the  radius  of  the  suction  tube,  is  smaller  than  one,  the  radius  of  the  vortex  core  is 
approximately  fixed  by  the  radius  of  the  vortex  tube.  When  the  diameter  of  the 
orifice  tube  is  small,  a  sudden  increase  of  the  vortex  diameter  occurs,  in  agreement 
Maxworthy's  observations  (1972).  The  latter  related  this  phenomena  to  vortex 
breakdown.  The  present  experimental  investigation  determines  the  range  of 
parameters  for  which  this  phenomenon  occurs. 


A  Hybrid  Vortex  Method  with  Deterministic  Diffusion 


Habib  N.  Najm 

Texas  Instruments,  SPDC,  Dallas,  TX 


In  the  general  class  of  particle  methods,  the  vortex  method  is  most  convenient  for  mod¬ 
eling  high  Reynolds  (Re)  number  vortex  flows.  This  is  particularly  due  to  the  following: 
(1)  the  Lagrangian  solution  of  the  Navier-Stokes  equations  eliminates  the  need  to  discretize 
the  non-linear  inertia  terms,  leading  to  good  numerical  stability  at  high  Re,  and  (2)  the 
restriction  of  computational  elements  to  the  regions  of  the  flow  exhibiting  shear  and  finite 
vorticity  leads  to  significant  numerical  efficiency. 

The  present  work  maintains  the  above  advantages  of  the  vortex  method  while  proposing 
an  accurate  scheme  for  modeling  the  diffusion  component  of  the  flow  equations.  The  method 
is  based  on  a  fractional  step  solution  of  the  two-dimensional  vorticity  transport  equation. 
The  advection  step  is  implemented  using  lagrangian  propagation  of  vortex  elements,  while 
the  diffusion  step  uses  a  hybrid  lagrangian-eulerian  implementation.  In  this  diffusion  scheme, 
both  the  vorticity  field,  u>(x,  y),  and  its  Laplacian,  V2o;(x,  y),  are  evaluated  from  the  summa¬ 
tion  of  the  individual  analytical  elemental  fields  at  any  point  in  the  computational  domain. 
In  particular,  both  quantities  are  computed  analytically  at  the  cell  centers  of  a  suitable  grid 
overlaid  on  the  computational  domain.  A  second  order  time  integration  of  the  diffusion 
equation  is  then  utilized  to  arrive  at  the  new  values  of  the  vorticity  at  the  cell  centers.  These 
are  used  to  generate  a  new  set  of  vortex  elements  at  each  time  step,  in  the  same  fashion 
utilized  in  the  initial  discretization  of  the  vorticity  field  at  time  t  =  0.  This  process  of 
continuous  regeneration  of  vortex  elements  at  each  time  step,  based  on  the  diffusion  of  the 
overall  vorticity  field,  allows  for  creation  of  new  elements  where  necessary  to  maintain  a  time 
accurate  representation  of  the  flow.  Further,  the  absence  of  grid  based  discretization  of  V2u; 
maintains  the  desireable  grid-free  solution  of  the  Navier-Stokes  equation  at  each  time  step. 

The  scheme  is  presented  and  demonstrated  on  an  incompressible  flow  problem  that 
involves  the  viscous  decay  of  a  vortex  tube  in  an  unbounded  two-dimensional  domain.  The 
analytical  solution  of  this  problem  is  known,  and  is  used  to  compute  the  error  associated  with 
the  numerical  solution,  to  study  its  convergence  characteristics,  and  to  compare  the  present 
results  against  those  available  from  the  random  vortex  method1,  and  the  deterministic  vortex 
method  due  to  Cottet  and  Gallic2,3.  Results  demonstrate  the  improvement  in  accuracy 
provided  by  this  hybrid  scheme  over  a  wide  range  of  Reynolds  number. 
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AXISYMMETRIC  VORTEX  SHEET  ROLL-UP 

Monika  Nitsche 
University  of  Michigan 

The  evolution  of  3-d  axisymmetric  vortex  sheets  is  studied  nu¬ 
merically  for  the  special  case  of  flow  without  swirl.  The  axisymmetric 
sheets  roll  up  at  their  boundaries  to  form  a  circular  vortex  ring.  The 
computations  use  a  vortex  blob  method  applied  to  axisymmetric  vor¬ 
tex  sheet  motion. 

Two  examples  are  studied.  The  first  one  is  the  evolution  of  an 
initially  flat  circular  vortex  sheet,  which  has  been  produced  by  setting 
a  circular  plate  impulsively  into  motion  in  an  ideal  fluid  and  then 
removing  the  plate.  The  vortex  sheet  rolls  up  into  a  ring.  We  study 
the  effect  of  the  numerical  parameters  upon  the  ring’s  propagation 
velocity,  radius,  and  vorticity  distribution.  The  computations  are 
compared  to  analytical  results  by  Taylor  (J.  Appl.  Phys,  vol.  24,  104) 
and  Saffman  (Stud,  in  Appl.  Math.,  Vol  XLIX,  No  4). 

The  second  example  is  the  formation  and  roll-up  of  a  vortex  sheet 
at  the  edge  of  a  circular  tube.  Well-documented  experiments  have 
been  performed  by  Didden  (ZAMP,  vol  30,  101),  in  which  a  piston 
drives  fluid  out  of  a  tube,  and  dye  particles  are  used  to  visualize 
the  vortex  sheet  formed.  We  compare  the  numerical  results  to  the 
experiments,  and  find  good  agreement  in  the  shape  of  the  spiral  and 
several  parameters  describing  it,  as  well  as  in  the  behaviour  in  time 
of  the  center  of  the  spiral. 


3-D  SIMULATION  OF  FREE  VORTEX  RINGS  AND  VORTEX 
RINGS  INTERACTING  WITH  A  WALL 


P.  Orlandi  and  R.  Verzicco 
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The  vortex  ring  is  a  basic  coherent  structure  and  the  detailed  comprehension  of  its  3- 
D  evolution  is  a  basic  feature  to  understand,  when  it  evolves  in  an  unbounded  domain,  the 
entrainment  of  external  fluid  into  round  turbulent  jets  and,  when  it  interacts  with  no-slip 
boundaries,  the  bursting  event  in  turbulent  boundary  layers.  The  aim  of  the  paper  is  to 
perform  finite  difference  numerical  simulations  second  order  accurate  in  time  and  in  space  to 
have  a  full  description  of  the  complex  vorticity  field. 

An  axisymmetric  calculation  has  been,  initially,  performed  in  both  cases  with  the  initial 
azimuthal  vorticity  u  =  -L_e-r2/ff3-  Because  the  Gaussian  distribution  is  not  a  solution  of  Euler 
equation,  in  the  case  of  a  free  ring,  the  study  has  been  focused  on  analysing  which  u/r  =  /(r/>) 
relationship,  (t/>  is  the  Stokes  streamfunction  in  a  coordinate  system  translating  with  the  ring) 
is  obtained  by  a  direct  numerical  simulation  performed  at  a  Re r  =  T/v  =  5500.  As  a  first  check 
of  the  precision  of  the  numerical  simulation  we  have  shown  that  the  calculated  dimensionless 
translation  velocity  of  the  ring  with  o  =  0.413,  Vt  =  2.23,  agrees  with  the  analytical  expression 
Vt  s=  2.49.  Moreover  vorticity  contour  plots  shows  that  immediately  after  the  ring  starts  its 
translatory  motion,  a  wake  is  generated;  this  wake  was  experimentally  observed  by  Maxworthy; 
from  the  simulation  we  conjectured  that  the  wake  is  formed  by  the  readjustment  of  the  initial 
vorticity  distribution  to  that  of  equilibrium,  for  which  scatter  plots  of  ui/r  =  fty)  show  a  good 
collapse  of  the  data.  We  obtained  a  function  /  which  smoothly  goes  from  0,  in  the  region 
external  to  the  core,  to  a  linear  function  in  the  region  within  the  core  of  the  ring. 

In  the  case  of  the  ring  interacting  with  a  non-slip  wall  we  performed  the  simulation  at  dif¬ 
ferent  Reynolds  numbers  and  we  compared  the  trajectories  of  the  centres  of  primary,  secondary 
and  tertiary  vortices  with  in  the  experiment  of  Walker  Tt  al.  Even  if  at  all  Reynolds  number, 
the  secondary  ring,  generated  at  the  wall,  causes  the  rebound  of  the  primary  ring  from  the  wall, 
only  at  high  Reynolds  number  a  pairing  between  secondary  and  tertiary  rings  forms  a  new  ring 
with  enough  circulation  to  migrate  far  from  the  wall  region. 

The  three-dimensional  simulation  of  the  Navier- Stokes  equations  for  the  free  ring,  with 
periodic  conditions  in  the  direction  of  translation,  has  been  performed  with  several  initial  per¬ 
turbations  in  the  azimuthal  direction.  In  all  cases  the  translation  velocity  does  not  change  from 
that  in  the  axisymmetric  case.  The  linear  stability  analysis  predicts  that,  for  this  translation 
velocity,  the  most  unstable  mode  is  that  for  n  =  5.  The  numerical  simulation  shows  that, 
giving  a  n  =  5  sine  waves  or  a  random  disturbance  with  the  same  amplitude  e/r  =  0.02  the 


perturbation  grows  in  time  maintaining  in  the  first  case  the  n  =  5  wave  (Fig. la)  and  in  the 
second  case  choosing  an  =  6  wave  (Fig.  lb).  A  further  case  has  been  performed  with  n  =  3  and 
we  observed  that  the  the  n  =  3  mode  is  dissipated  and  the  n  =  6  mode  is  that  which  grows. 

By  the  direct  simulation  we  proved  that  for  thick  free  rings  the  fundamental  and  higher 
modes  are  excited  but  only  the  most  unstable  persists  and  grows  in  time.  In  all  cases  flow 
visualisation  of  isolevel  surfaces  of  u>e  and  u>a  =  show  the  deformation  of  the  ring 

structure.  u2  gives  the  greater  contribution  to  uia  and  uz  initially  is  located  in  very  long  and 
thin  structures  at  the  interior  of  the  ring.  These  patches  of  opposite  sign  vorticity  are  generated 
near  the  centre  of  the  ring  core,  are  deformed  by  the  wake  and  become  elongated  structures 
which  by  mutual  induction  migrate  towards  the  centre  where  are  dissipated. 

In  the  3-D  simulation  of  the  vortex  ring  impacting  a  solid  wall  we  observed  that  azimuthal 
instabilities  are  effecting  mainly  the  secondary  ring,  which  is  advected  upwards  by  the  primary 
ring;  in  this  first  stage  (Fig.2a)  the  initial  n  =  5  perturbation  grows.  In  a  second  stage,  while 
the  secondary  ring  is  advected  and  compressed  at  the  interior  of  the  primary  ring  (Fig. 2b),  the 
radius  decreases  and  a  well  defined  n  —  10  mode  is  observed.  When  the  secondary  ring  is  very 
close  to  the  wall  a  substantial  amount  of  u>z  is  generated  and  five  rings  of  smaller  dimensions 
are  generated  that  interact  and  pair  with  the  tertiary  ring  during  its  translation  towards  the 
centre  (Fig.2c).  The  tertiary  ring  and  the  smaller  new  rings  pair  and  produce  a  structure  with 
enough  circulation  to  move  far  from  the  wall  (Fig.2d).  This  complex  interaction  of  secondary 
vorticity  patches  induces  on  the  primary  ring  a  n  =  10  deformation.  For  the  case  of  vortex 
rings  interacting  with  solid  wall  there  was  no  any  theoretical  analysis  but  we  found  that  our 
vorticity  magnitude  contour  plots  resemble  very  close  the  flow  visualisation  of  Walker  et  al.  At 
the  workshop  comparison  between  the  ejection  velocity  in  the  axisymmetric  and  that  in  the 
3-D  case  will  be  presented,  together  with  the  differences  when  a  random  initial  perturbation  is 
given.  Furthermore  will  be  presented  vorticity,  and  r.m.s.  velocity  distributions  which  are  very 
difficult  to  measure. 
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Abstract 


Vortex  dipoles  are  common  features  of  geophysical  flows,  and  they  play  an  important  role 
in  their  dynamics.  In  the  ocean,  dipolar  vortices  may  be  generated  in  various  ways  (e.g.  as  a 
results  of  shedding  from  unstable  boundary  currents  or  due  to  localized  wind  forcing)  and  they 
provide  an  important  mechanism  in  the  transport  of  physical  properties  (heat,  pollutants  ).  In 
the  atmosphere,  dipolar  flow  structures  in  the  form  of  blocking  systems  tend  to  have  a  stabilizing 
influence  on  the  local  weather.  Within  the  context  of  stability  of  such  flow  structures  it  is  of 
importance  to  know  whether  the  structure,  once  perturbed,  relaxes  towards  its  initial  (stable) 
state.  It  is  easy  to  show  that  any  functional  relationship  u>  =  f(ip)  between  the  vorticity  u;  and 
the  streamfunction  ip  satisfies  J(uj,  ip)  =  0  and  is  thus  a  stationary  solution  of  the  inviscid  Euler 
equations.  In  the  present  study  we  consider  the  Lamb  dipole,  which  travels  with  a  constant 
velocity  while  preserving  its  shape.  For  this  dipolar  vortex  the  relation  between  uj  and  ip,  in  a 
reference  frame  travelling  with  the  dipole,  is  a  linear  one.  i.e.  u>  =  k2ip. 

Analytical  approaches  based  on  linear  approximations  are  not  very  helpful  to  study  the 
stability  of  dipoles  as  has  been  recently  shown  by  Nycander  (1992)  for  modons,  which  are 
dipoles  conserving  their  properties  in  a  beta  plane.  Numerical  approaches  are  more  useful,  in 
fact  different  types  of  perturbations  can  be  applied  to  see  in  which  case  the  dipole  reaches  a 
new  steady  configuration  and  which  functional  relationship  f(ip)  is  yielded.  We  have  performed 
numerical  simulations  with  two  different  initial  configurations  out  of  equilibrium.  In  the  first 
one  the  vorticity  is  confined  within  an  ellipse  with  the  longer  axis  aligned  with  the  direction  of 
translation  and  in  the  second  one  the  longer  axis  is  orthogonal  to  the  translational  direction.  In 
both  cases  we  observed  that,  at  Re  =  10000,  a  value  very  close  to  the  inviscid  conditions,  the 
structure  in  a  first  stage  sheds  a  tail  and  later  on  forms  a  steady  dipole.  The  time  necessary  to 
reach  the  steady  configuration  depends  on  the  disturbance.  In  the  case  of  an  ellipse  elongated 
in  the  direction  of  translation  the  same  initial  f(ip)  is  obtained,  (Fig. la)  while  in  the  other  case 
a  relationship  with  two  slopes  is  reached  (Fig. lb).  The  reason  of  the  reduced  slope  is  due  to 
the  greater  distance  between  the  two  vorticity  peaks  than  those  for  the  Lamb  dipole.  Couder 
and  Basdevant  (1986)  found  a  similar  flattening  near  the  axis  in  their  numerical  simulations. 
The  core  of  the  dipole  on  the  contrary  is  well  restored. 

A  different  perturbation  has  been  given  by  letting  two  dipoles  to  collide  at  different  angles, 
as  done  in  the  experiment  by  van  Heijst  and  Flor  (1989).  In  this  case  during  the  impact  the 


dipoles  are  deformed,  interchange  partner  and  finally  travel  again  along  a  straight  trajectory. 
The  numerical  simulation  has  been  performed  by  assuming  that  the  dipole  is  impacting  with  a 
free-slip  wall.  We  considered  the  collision  at  three  different  angles  90,60  and  30  degrees.  The 
scatter  plots  are  given  in  Fig.2  a-c.  Fig.2a  shows  that  for  the  collision  at  90°  the  law  f{ip)  —  k  ip 
is  reached.  On  the  other  hand  Figs.2b  and  2c  show  that,  when  the  dipole  impacts  the  wall  at 
an  angle  a  functional  relationship  with  a  discontinuity  on  the  slope  is  obtained  and  that  the 
region  responsible  of  the  higher  gradient  is  that  close  to  the  wall.  To  obtain  these  plots  we 
analysed  the  field  of  the  dipole  which  was  closer  to  the  wall  and  was  moving  faster. 

From  these  two  different  perturbations  we  conclude  that  the  region  close  to  the  peak  vor- 
ticity  is  very  stable,  while  the  region  close  to  the  centre  of  the  dipole  varies  depending  on  the 
perturbation  given.  Dipoles  with  a  flat  region  are  those  that  can  easily  loose  the  coherence  and 
split  into  two  monopoles. 

We  were  also  interested  to  understand  which  is  the  global  quantity  representing,  in  the 
best  way,  when  the  steady  state  is  reached.  We  have  thus  calculated,  in  the  region  with  the 
dipolar  structure,  the  excess  energy,  the  enstrophy  and  the  integral  of  ulog(u>)  that  represents 
the  entrophy  of  the  flow  field.  We  found  that  when  the  minimum  for  each  quantities  is  reached 
a  well  defined  f(ip)  is  observed.  Time  histories  of  these  quantities  will  be  presented  at  the 
conference. 
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MODELISATION  D'UN  TOURBILLON 
EN  ECOULEHENT  TURBULENT  INCOMPRESSIBLE 


par  Philippe  PASCAL 
Resume 

Les  icoulements  resultant  de  Involution  d'une  nappe  fluide,  engendrie 
par  un  dicollement,  donnent  naissance  a  des  structures  tourbillonnaires 
enroulies,  que  1'on  nomme  d'une  man i ire  plus  courante:  tourbi lions. 

En  Airodynamique,  ces  tourbillons  peuvent  avoir  des  repercussions 
importantes  sur  les  conditions  de  vol  a  grande  incidence  des  avions  et  des 
engins  munis  d'ailes  4  forte  fleche.  En  effet,  sous  1 'action  des  champs  de 
pression  associes  a  la  voilure,  ces  tourbillons  sont  susceptibles  de  subir  une 
disorganisation  brutale,  plus  connue  sous  le  nom  d'eclatement  tourbi 11 onnai re. 
Ce  phinomine  entratne  une  modification  profonde  de  1 'icoulement  sur  l'aile, 
dont  la  provision  reste  tres  incertaine.  L'objet  de  ce  travail  est  de  proposer 
une  mithode  de  calcul  capable  de  deer ire  ce  type  d' icoulement.  Une  itude 

thiorique  a  done  iti  rialisie,  dans  le  but  de  comparer  les  diffirentes 

mithodes  de  priditermi nation  de  1'ivolution  d'un  tourbillon,  dans  son 

processus  d'eclatement.  Les  mithodes  proposies  sont  basies  sur  l'hypothese 

d'axisymitrie  de  1 'icoulement,  et  sur  un  traitement  statistique  classique  de 
la  turbulence  avec  fermeture  en  un  point.  Trois  approches  diffirentes  ont  iti 
menies: 

-  en  utilisant  les  iquations  statistiques  du  mouvement  icrites  avec  les 
approximations  de  la  couche  limite,  et  en  introduisant  un  model e  de 
turbulence  algibrique, 

-  en  appliquant  les  mimes  iquations,  mais  en  utilisant  l'A.S.M, 

-  en  risolvant  directement  les  iquations  statistiques,  sans  approximation  et 
avec  un  modile  algibrique.  Une  confrontation  des  risultats  du  calcul  avec 
certaines  donnies  expirimentales  permet  de  virifier  1'adiquation  entre  le 
calcul  et  1 'expirience,  ainsi  que  les  tendances  sur  la  validiti  des  modiles  de 
turbulence. 

Une  itude  expirimentale  a  iti  rialisie  pour  mieux  apprihender  les 
diffirents  phinomines  constituent  un  tourbillon,  et  notament  analyser  et 
quantifier  1' influence  de  la  turbulence  sur  le  diveloppement  d'un  tourbillon 
stable  non  iclati.  Pour  cela,  on  a  utilisi  un  dispositif  permettant  de 
produire  un  tourbillon  de  bout  d'aile.  La  structure  interne  du  tourbillon 
prisente,  au  centre,  une  survitesse  qui  tend  k  diminuer  lorsque  le  tourbillon 
ivolue  dans  sa  direction  axiale.  Le  champ  de  vitesse  extirieur  est  bien 
reprisenti  par  un  tourbillon  potentiel.  Toutefois,  on  a  remarqui  que  les 
mouvements  latiraux  d'ensemble  de  1 'icoulement  doivent  etre  pris  en  compte 
dans  1 'interpritation  des  risultats  expirimentaux,  car  ils  induisent  des 
fluctuations  k  basse  friquence,  qui  se  superposent  k  la  turbulence  proprement 
dite. 


A  slightly  diffusive  Contour  Dynamics 
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The  use  of  a  piecewise  constant  approximation  of  the  vorticity  field  gives  a  simple  and  fruitful  approach 
to  the  Lagrangian  simulation  of  the  coherent  structures  motion  in  a  two-dimensional  intnscid  context.  The 
way  is  a  straightforward  generalization  of  the  Zabusky  idea  to  multi-levels  distribution  of  vorticity  based  on 
the  linearity  of  the  vorticity-velocity  rule.  In  the  following  we  outline  a  natural  extension  of  the  Contour 
Dynamics  to  a  slightly  diffusive  fluid  by  using  a  simple  approximation  of  the  vorticity  fluxes  between  adjacent 
contours. 

The  velocity  field  induced  by  a  piecewise  constant  vorticity 

R 

ut  =  53H.XIMO 
»=1 

is  given  in  terms  of  the  contours  {dDj}j=\ . „  and  vorticity  jumps  {Hj}y=i . n  across  these  by 

««(*)  =  ~  f  . .  t{s)G(x  -  y)ds 


in  which  G{x )  =  — log\x\.  An  approximate  description  of  the  vorticity  dynamics  is  given  by  using  a  straight 

panel  discretization  of  each  contour  dDj  with  a  suitable  panel  duplication  procedure. 

The  fascinating  aspect  of  this  method  is  the  absence  of  any  form  of  numerical  viscosity  that  enables  to 
follow  the  small  scales  motion  as  clearly  shown  in  the  collision  of  Lamb  dipoles  of  Fig.  1.  This  suggests  to 
modify  the  Contour  Dynamics  technique  in  order  to  study  the  small  scales  motion  in  presence  of  a  slight 
diffusion  due  to  real  viscosity.  A  simple  way  is  to  consider  the  vorticity  jumps  as  time  dependent 
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(we  Msnm<>  Do  D  D\  D  •  •  •  D  Dn)  and  to  integrate  the  Helmholtz  equation  in  all  the  annular  regions 
Di- 1  —  Di,  i  —  1, . . . ,  n  and  also  in  D„.  By  using  a  suitable  approximation  of  the  vorticity  fluxes  exchanged 

(  d[w]i  1 

by  adjacent  levels  finally  we  obtain  a  first  order  differential  system  in  terms  of  <  —  —  -  > 

1  “*  J  »=!,.. .,n 


The  accuracy  of  the  time  integration  as  well  as  the  effects  on  the  numerical  solution  of  the  vorticity 
fluxes  approximation  are  evaluated  computing  some  “global  quantity”  of  the  flow  that  in  an  inviscid  context 
play  the  role  of  first  integrals  of  the  motion.  These  computations  can  be  carried  out  in  a  simple  way  as  long 
as  the  limit  T  of  the  circulation  around  spheres  Br  for  R  — *  oo  is  numerically  retained.  The  most  useful 

relations  involve  the  second  moment  of  vorticity  I  =  \x\ 2u>dx 

Jr * 

m  =  m + ^ri 


and  the  excess  energy  £r[t)  —  ~  f  tpudx  for  which  we  have 

2  Jri 

id') = f/(°)  -  i 


where  E  is  the  enstrophy  of  the  flow. 

As  an  example  of  calculation  of  the  previous  quantities  we  show  in  Fig. 2  the  axisymmetrization  of  an 
elliptical  non  uniform  vortex  in  an  inviscid  and  in  a  viscous  flow.  In  this  case  the  effects  of  the  vorticity 
fluxes  approximation  appear  clearly  negligible. 


Estimation  numerique  d’ecoulements 
tridimensionnels  avec  une  methode  particulaire  a 

poids  constant 

E.  Rivoalen 


1  Introduction 

L’objectif  de  c.e  travail  est  la  mi.se  au  point  de  methodes  niuneriques  susceptibles  detre 
utilisees  pour  calculer  revolution  d’un  tourbillon  emis  d’un  profil  portant  sur  des  dis¬ 
tances  importantes.  En  effet  le  calcul  de  transport  d’un  sillage  tourbillonnaire  pendant 
des  temps  relativeinent  grands  pose  des  problemes  a  la  fois  de  precision  des  algorit limes 
utilises  et  de  viscosite  numerique.  Les  imprecisions  peuvent  etre  acceptables  si  Ton  ne 
s’interesse  qu’aux  caracteristiques  hydraudynamiques  du  profil  generateur  de  sillage  , 
mais  ne  le  sont  plus  si  l’on  regarde  l’interaction  de  ce  sillage  avec  d’autres  corps  comme 
par  exemple  une  helice,  ime  carene  de  bateau  ou  encore  ime  surface  libre.  Quand  on 
s’interesse  au  transport  du  sillage  sur  des  distances  relativeinent  grandes  il  apparait 
egalement  necessaire  de  prendre  en  compte  les  phenomenes  de  diffusion  si  I’on  veut 
doimer  une  description  realiste  de  l’ecoulement.  Nous  presentons  dans  cette  etude 
une  moderation  de  l’ecoulement  tourbillonnaire  tridimensiomielle  par  la  methode 
des  tourbillons  ponctuels.  La  necessity  de  conserver  le  tourbillon  sur  de  tres  longues 
distances  nous  amene  a  travailler  avec  une  methode  utilisant  une  description  lagrang- 
ienne  de  l’ecoulement.  Nous  utilisons  ici  le  cadre  de  la  methode  de  discretisation 
particulaire  tridimensiomielle  mise  au  point  par  Relibach  en  1977. 

2  Transport  du  tourbillon  a  poids  constant 

L’objectif  numerique  est  de  construire  une  methode  d’integration  de  l’equation  de 
Helmholtz.  Dans  ce  cas  1’equation  de  transport  du  toiu-billon  comporte  un  terme 
supplementaire,  le  terme  de  deformation.  Ce  dernier  est  decompose  en  un  champ 
tangent  au  tourbillon  local  qui  represente  l’elongation  du  tourbillon  et  un  champ  or¬ 
thogonal  qui  represente  une  reorientation  du  tourbillon  local  done  une  transformation 
a  module  constant.  On  remarque  que  le  terme  tangent  peut  sinterpreter  comme  une 
redistribution  des  particules  le  long  des  filets  tombillonnaires  du  fait  des  elongations 
et  contraction  que  l’ecoulement  externe  lui  fait  subir.  Ce  terme  apparaissait  deja  dans 
le  terme  eoveetif  de  l’equation  de  transport,  il  disparait  done  de  cette  equation.  Sous 
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forme  discrete  1’equation  de  transport  du  vecteur  tourbillon  s’ecrit  : 


1  §  = 

Cette  decomposition  revient  a  ne  considerer  localement  que  le  champ  de  vitesse  or¬ 
thogonal  au  tourbillon. 


3  Modelisation  de  la  diffusion 

La  prise  en  compte  des  phenomenes  de  diffusion  est  traitee  comme  un  reamenagement 
de  la  repartition  des  points  de  discretisation.  Nous  proposons  une  methode  mixte  dans 
laquelle  la  composante  tangente  du  vecteur  tourbillon  est  diffusee  par  une  methode 
de  depacement  detenniniste  (l’equation  de  diffusion  pure  est  couplee  a  mie  equation 
de  transport)  et  les  deux  autre  composantescomme  une  rotation  du  vecteur  tourbil¬ 
lon  (transformation  a  module  constant).  Ce  modele  permet  de  rendre  compte  de 
revolution  de  la  vortite  sous  l’effet  de  la  diffusion  et  en  particulier  de  l’expension  d’un 
sillage. 
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Un  modele  d'ecoulement  turbulent  bidimensionnel 


Raoul  Robert 

Universite  Claude  Bernard 
69000  Villeurbanne  Cedex 


On  utilise  la  theorie  des  etats  d'equilibre  statistiques  de  l'equation  d'Euler  2D  pour  construire  un 
modele  de  turbulence  en  evolution  pour  un  fluide  bidimensionnel  incompressible  faiblement 
visqueux. 


THE  USE  OF  THE  “NATOKAP”  STYLE  FILE  FOR 
DELIVERING  MANUSCRIPTS  IN  MTEX 


The  use  of  the  NATOKAP  style  file  for  delivering  manuscripts  in  IXIjrX: 

1.  Use  standard  LaTeX.  Anything  allowed  in  LaTeX  may  be  used. 

Begin  your  article  with: 

\documentstyle{IATOKAP} 

The  following  items  are  different: 

a.  ‘Opening-environment. 

This  environment  helps  to  print  the  title,  author,  affiliation  and  date  in_the  required 
format. 

Use: 


\begin  {opening} 

\title  {title  of  your  paper}  %  This  is  a  new  command 

\subtitle  {if  there  is  any}  !•  This  is  a  new  command 

\author  {one  author's  name} 

\institute  {address  of  your  institute}  '/,  This  is  a  new  command 
Vend  {opening} 

b.  Do  not  use  the  command. 

c.  Acknowledgements. 

The  acknowledgements  command  can  be  used  for  making  these.  The  text  can  be 
written  by  starting  on  the  next  line. 

d.  The  commands: 

Varticle 

\book 

Vinproceedings 

VbibitemO  7.  The  brackets  must  be  empty! 
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have  been  implemented  to  ease  the  use  of  bibliographical  data.  They  have  been 
‘borrowed’  from  BibTEX,  so  BibTEX  need  not  to  be  used.  If  the  standard  items 
required  by  the  first  three  commands  are  not  available  (see  your  LaTeX  manual), 
you  are  kindly  requested  to  use  bibitea  and  do  the  formatting  yourself  (e.g.  with 
titles  of  articles). 

N.B.  If  you  would  like  to  use  BibTEX  then  you  are  invited  to  do  so,  but  you  need 
to  indicate  where  the  BibTEX  files  are. 

If  you  encounter  problems,  please  call:  (78)  334  316 

or  send  a  message  to  email:  SURF406@KUB.NL 

Your  document  should  be  structured  along  the  following  lines: 

\documentstyle{natokap} 

V.  The  'OPEIIIG' -environment 
\begin{opening> 

\titleO 

\aut hor { Author 1 } 

\author{Author2} 

Ninstitute-Clnstitute  of  Authorl  and  Author2> 

\author{Author3} 

\inatitute{Institute  of  Author3> 

Y.  etcetera 
\end{opening} 

\begin{docuaent> 

\begin{ abstract} 

Tekst  of  abstract. 

\endfabstract} 

\section{The  Sectiontitle  Should  lot  Be  Typed  with  Capitals  Only} 
\begin{f igur e} 

\vspace{??cm}  V,  You  can  indicate  here  the  space  your  figure  needs 
NcaptionQ {caption  text} 

\end{f igure} 

Nbegin-Ctable} 

\capt ion [] {caption  text} 

%  Here  the  data  and  structure  of  the  table 
Vendftable} 

%  I.B.:  You  aay  use  the  "tabular"  environment  as  veil 
Nacknovledgements 

If  you  would  like  to  thank  people  or  organizations 
\beg  in{  t  hebibl  iographyH} 

%  Space  between  second  pair  of  brackets  is  empty 
Nendfthebibliography} 


Nendfdocument} 
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ABSTRACT.  These  instructions  are  intended  to  help  and  guide  contributing  authors  and  their  typists  when 
preparing  articles  as  camera-ready  output  on  a  laserprinter. 


1.  Format  Instructions 

1.1.  INSTRUCTIONS  FOR  TYPISTS 

Tins  document  has  been  output  on  a  QMS  PS-800  laserprinter  using  the  software  Flexicon  with 
justified  text,  automatic  centering,  italics,  and  boldface  facilities.  Margins  were  set  to  a  width  of 
15.2  cm  (36  picas),  and  each  page  was  typed  in  Times  1 1  points  letter  for  the  main  text  with  a  13 
points  base-to-base  spacing  of  the  lines.  The  small  typeface  was  typed  in  Times  10  points  letter 
with  a  1 1  points  base-to-base  spacing  of  the  lines.  Each  page  was  typed  to  a  page  depth  of  21.6 
cm  (51  picas,  47  text  lines).  The  headings  used  are: 

2.  First-Order  Headings 


These  arc  typed  in  bold,  upper  and  lowercase  letters,  with  two  lines  of  space  above  the  heading, 
and  one  line  below.  The  text  after  the  heading  begins  at  the  left-hand  margin  (i.e.  it  is  not  in¬ 
dented  as  for  new  paragraphs). 

2.1.  SECOND-ORDER  HEADINGS 

These  are  typed  in  light  small  capitals,  at  the  left-hand  margin.  There  is  one  line  of  space  above, 
and  one  line  of  space  below  the  heading. 

2.1.1.  Third-Order  Headings.  Typeset  in  italics,  with  capital  initial  letters.  Run  the  following  text 
on  in  the  same  line,  and  leave  one  blank  line  above. 

New  paragraphs  are  to  be  indented  1  pica  from  the  left-hand  margin,  with  no  extra  space 
between  the  paragraphs. 

Whatever  your  final  decision  on  the  layout  of  headings,  etc.,  please  ensure  that  you  type  each 
page  to  the  correct  dimensions,  as  found  on  this  model.  Set  the  word  processor  margins  to  the 
appropriate  width,  and  type  to  the  same  overall  depth  on  each  page. 


3.  Paper 


If  you  are  using  a  word  processor  to  produce  your  text  on  a  laserprinter,  you  can  use  the  normal 
A4  (210  x  297  mm)  standard  paper.  The  layout  should  be  according  to  this  model  and  the 
instructions  mentioned  in  Section  1 . 1 . 


4.  Corrections 

It  is  assumed  that,  if  you  are  using  a  word  processor,  you  will  re-output,  where  necessary,  after 
editing  and  correction  on-screen. 

5.  Special  Problems 

5.1.  SPECIAL  CHARACTERS  AND  SYMBOLS 

It  can  happen  that  some  equations  or  formulae  in  your  work  cannot  be  easily  typed  because  they 
require  the  use  of  symbols  which  are  not  available  on  the  apparatus  you  are  using.  Some  common 
equipment,  for  example,  will  not  handle  the  accents  necessary  for  typing  most  European 
languages.  A  decision  will  have  to  be  made,  therefore,  on  how  to  treat  these  occurrences.  Ob¬ 
viously,  this  is  a  matter  of  degree:  if  you  are  only  citing  names  and/or  titles  of  non-English- 
language  works,  for  instance,  it  is  possible  and  perfectly  acceptable  to  insert  accents  by  hand  on 
the  finished  typescript,  or  to  leave  them  out  altogether. 

For  most  mathematical  and/or  logical  writings,  there  often  exists  an  alternative  formulation 
which  avoids  the  use  of  ‘difficult’  characters.  In  cases  where  it  is  genuinely  unavoidable,  there  is 
a  wide  range  of  unusual  symbols  and  letters  available  on  dry-transfer  lettering.  Handwritten 
equations  are  strongly  discouraged. 


